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Relative methods, which are performed with the assistance of reference materials, are widely used in photon 
activation analysis (PAA). On the contrary, absolute methods, which are conducted without any reference mate- 
rial, are rarely applied due to the difficulty in obtaining photon flux. To realize absolute measurement in PAA, 
we retrieve photon flux in the sample via Monte Carlo simulation and raise a novel procedure—quasi-absolute 
method. With simulated photon flux and cross section data from existing databases, it is possible to calculate 
the concentration of target elements in the sample straightforwardly. A controlled experiment indicates that re- 
sults from the quasi-absolute method for certain elements are nearly comparable to relative methods in practice. 
This technique of absolute measurement has room for improvement in the future and can serve as a validation 


technique for experimental data on cross sections as well. 
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I. INTRODUCTION 


Photon activation analysis (PAA), in most cases nowadays, 
is an accelerator-based radioanalytical technique. After the ir- 
radiation of the samples with high energy photons created by 
bremsstrahlung radiation of electrons, qualitative and quanti- 
tative information on the target elements can be obtained by 
detection of nuclear emissions following the irradiation with 
LINAC. The first article of PAA dated back to 1951 when 
Gaudin and Pannell of MIT tried to determine the amount of 
beryllium in low grade beryl ores by photodisintegration of 
Beryllium [1]. They commented on this method as “rapid, 
simple, and nondestructive”. In 1954, Basile proposed to 
analyze some of the light elements using photonuclear reac- 
tion induced by bremsstrahlung radiation from a betatron [2]. 
His paper led to the wide use of the bremsstrahlung radia- 
tion as the high energy photon source of PAA. A chapter in 
Encyclopedia of Analytical Chemistry gave a brief but com- 
plete summary of PAA research [3]. In recent years, PAA 
has expanded its applications in radiotherapy, meteorology, 
geochemistry, archeology, industrial material, environmental 
studies, etc. [4-10]. 

Figure | shows the basic principle of PAA. Photon acti- 
vation usually occurs in the giant dipole resonance (GDR) 
region of the target nucleus. A high energy photon above 
the particle emission threshold enters the nucleus and inter- 
acts with nucleons. Then the excited nucleus equilibrates and 
turns into a statistical “compound nucleus” state. This “com- 
pound nucleus” will de-excite through residual strong inter- 
action by emitting prompt gamma rays, neutrons, protons, 
and even alpha particles. Following the “compound” nucleus 
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stage, the product nucleus generally is still unstable because it 
is usually a proton-rich (or neutron-rich) nucleus. It continues 
to de-excite to the final nucleus through weak interactions, 
such as electron capture and beta decay. Delayed gamma 
rays are emitted after those weak interactions. In most cases 
of PAA study, researchers concentrated on delayed gamma 
ray emission with a relatively long half-life. Those gamma 
rays can be detected by HpGe spectrometers and generated 
spectra with the assistance from computers. Hence, the pro- 
cedure of PAA is quite straightforward: from the spectra of 
delayed gamma rays, one acquires their energy lines and in- 
tensity, which corresponds to certain product nuclides. From 
the information on the products nuclides, one can trace back 
to quantity and species of the target nuclides according to the 
corresponding nuclear reaction channels. 


IL CALCULATIONS 


According to Ref. [11-13], the net peak area (or counts) P 
of the characteristic gamma line in the spectrum is 


P =u | A(t)dt 


mhL 
“ee (1 _ ei) (1 = e rte e—Ata (1) 
Emar 
x | y(E)o(E)dE, 
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where 77 is the detector efficiency, 0 is the branching ratio of 
the reaction channel, ¢ is the absolute intensity of the gamma 
line, À is the activity of radioactive nuclide, m is the mass 
of the total sample, Cm is the concentration of the nuclides of 
interest, h is the natural abundance of the target nuclide, L 
is the Avogadro constant, À is the decay constant of product 
nuclide, À, is the atomic mass of the target nuclide, y(£) is 


050201-1 


202306.00459v1 


chinaXiv 


SUN Zai-Jing et al. 


neutron,proton, alpha 


l l 
ee6 Y 


ChinaXiv&@E#A 


-H 


Nucl. Sci. Tech. 25, 050201 (2014) 


delayed 
gamma 


pe Oe BU 


target nucleus “compoud nucleus” 


product nucleus 


final nucleus 


Fig. 1. (Color online) Basic principle of PAA. 


the energy differential photon flux, o(F) is the cross section 

of the photonuclear reaction of interest, t; is the irradiation 

time, te is the counting time, and tą is the decay time from 

the end of irradiation to the start of spectrum collection. 
Equation (1) can be rearranged as 


E PA,X 
Cm nOCMAL(A — es) (1 — ete je Ma 
i (2) 


(E)o(E)àE 
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On the left side of Eq. (2) is the concentration €m of a cer- 
tain element in the sample, which is the ultimate value PAA 
seeks. On the right side mostly are the parameters we can 
measure (P, 7, M, ti, te, ta), or parameters we can find the 
values in existing databases (0, ¢, h, L, À, Ar, o(E)). The 
only parameter we normally cannot measure directly is the 
energy differential photon flux y(£). 

How to find the energy differential photon flux y(F)? As 
mentioned before, we need assistance from computer simula- 
tions. From the simulated photon yield, it is feasible to derive 
photon flux and then calculate the concentration of the target 
nuclide directly based on the Eq. (2). 

The output of the computer simulation is the photon yield 
Y(£). To obtain photon flux y(£), one needs to multiply it 
with average beam current, Jjeam, Of the LINAC as below. 

PIE) = Y(E) x =, 6) 
where qe is the charge of an electron. Inserting Eq. (3) into 
Eq. (2), one gets 


7 PA,AGe 
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To simplify Eq. (4), we define the integral in the denomi- 

nator as the reduced reaction rate Ip, the Eq. (4) evolves into 
E PArAqe 

Cm = ndCmhL(1 = ei) (1 = e— Ate je Ald ThoamT R 


(5) 


In Eq. (5), the uncertainty from atomic mass A,., decay 
constant A, electron charge qe, detector efficiency 7, branch- 
ing ratio 0, absolute intensity Ç, natural abundance h, Avo- 
gadro’s number L, and time parameters (t;, te, ta) usually 
can be ignored. The primary source of uncertainty comes 
from uncertainty in the net peak area P, mass of sample m, 
beam current Jbeam, and reduced reaction rate IR. Hence, the 
uncertainty propagation of Eq. (5) follows 


2 2 2 2 
TR (+) (22) (<=) (=) © 
P m Theam IR 
One should note that the uncertainty of the reduced reac- 


tion rate in Eq. (6) is the combination of the uncertainties of 
photon yield and cross sections. 


IHI. EXPERIMENTAL SETUP, SIMULATIONS, CROSS 
SECITONS, AND PHOTON YIELD 


A. Experimental setup 


The experiment to realize quasi-absolute method was done 
by the 44MeV L-band pulsed LINAC at the Idaho Accel- 
erator Center. In the experiment, the peak current is about 
300 mA, pulse width is 2.2 us, repetition rate is 120 Hz, and 
the peak energy is about 30 MeV. The side view of the exper- 
imental setup is shown in Fig. 2. 

Electrons are first produced by the hot cathode and then 
accelerated by a series of alternating RF electric fields in the 
acceleration cells. They form a focused electron beam from 
magnetic fields. The radius of the electron beam is around 3 
millimeters. The energy distribution of the electrons (Fig. 3) 
was measured with a Faraday cup, a magnetic spectrome- 
ter and beam-split to allow quasi-monochromatic electrons to 
pass after magnetic fields bend the track of electrons. 

After the electron beam passes through a water cooling sys- 
tem and a path of air, it is absorbed by a 3 mm-thick tantalum 
electron-photon converter (or radiator). Bremsstrahlung radi- 
ation is the primarily physical process in the radiator. After 
the radiator, the electron beam almost completely stops and 
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Fig. 2. (Color online) Experimental setup of quasi-absolute method. 
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Fig. 3. (Color online) Energy distribution of the electron beam of 
LINAC at 30 MeV peak energy. 


the newly generated photon beam strikes an aluminum hard- 
ener of about 3 inches thickness. The hardener absorbs the 
residual electrons and ensures that the beam out of the radia- 
tor is entirely made of high energy photons. 

The urban particulate matter —standard reference material 
1648a [14]—from NIST was chosen as the sample for veri- 
fying quasi-absolute method. It was selected because it is a 
well-known certified reference material and widely used in 
instrumental analytical methods. The sample was wrapped 
into a 1 cm by 1 cm square area with Aluminum foil and put 
behind the hardener along the beam axis. The thickness of 
the sample is about 1 mm. 


B. Photon beam simulations 


Photon flux simulation was exactly based on the geom- 
etry and material of the experimental setup in Fig. 2. The 
simulated photon flux is on the surface of the sample. Sam- 
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Fig. 4. (Color online) Geant4 simulation of “photon shower” in ex- 
periments for quasi-absolute method (view point 0=45°, y=45°). 
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Fig. 5. (Color online) Comparison of the simulated electron-photon 
yield with Geant4 and MCNPX (log scale). 


ple matrix effects (e.g., in sample Compton scattering, y in- 
duced further nuclear reactions, etc.) are not considered since 
the sample is small and thin. Geant4 is used as our primar- 
ily simulation software. MCNPX applied as an accessorial 
validation. The Geant4 toolkit we applied to the simulation 
is version 4.9.6 [15]. The class library of the high energy 
physics (CLHEP) used in this version is 2.1.3.1 [16]. The 
Geant4 toolkit is installed on a computer server running 64 
bits GNU/Linux. The physics list of simulations includes all 
the electromagnetic processes, such as ionization, multiple 
scattering, positron annihilation, photoelectric effects, Comp- 
ton scattering, pair production, and bremsstrahlung radiation. 
In the simulations, the total number of electrons we applied 
was 10 million. The statistical uncertainty from the simula- 
tion was small since we processed such a large number of 
electrons. 
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Figure 4 depicts the Geant4 simulation of a photon beam 
for the quasi-absolute method in PAA. In the figure, the 
yellow dot is the tracking point. The relationship of the 
track color and its corresponding particle is: photon-green, 
electron-red, positron—blue, neutron-yellow. The yellow 
square disk at the end of beam line is the sample. 

Figure 5 shows the difference between simulated photon 
yields from Geant4 and those from MCNPX with the same 
geometry and materials shown in Fig. 2. The two simulated 
results are consistent but with some slight variance. The per- 
centage of differences is approximately 2% for the whole en- 
ergy range (0-30 MeV) and about 7.6% for the energy range 
of interest (10-30 MeV) of PAA. The source of this discrep- 
ancy may originate from the different cross section libraries 
and the different algorithms used to implement calculation of 
electromagnetic processes in the two simulation codes. 


C. Selection of cross sections 


Ten photonuclear reactions were selected for realizing the 
quasi-absolute method (see Table 2). They were chosen based 
on the following criteria: (1) concentrations of their target 
nuclides are certified; (2) the target nuclides in this sam- 
ple have a wide atomic number distribution: from light el- 
ements to heavy ones; and (3) their product nuclides have 
clear interference-free energy lines (pile-up can be ignored at 
the low counting rates in these experiments); (4) most of them 
are (y, n) reactions. For the same isotope, the (y, n) reaction 
channel has the largest cross sections in GDR region. 

Original data of cross sections is in exchange format (EX- 
FOR), which contains an extensive compilation of experi- 
mental nuclear reaction data [17-37]. These records were 
retrieved and validated via a consortium of organizations, in- 
cluding the International Atomic Energy Agency (IAEA), the 
National Nuclear Data Center (NNDC), and the Center for 
Photonuclear Experiments (CDFE). Our criteria for choos- 
ing cross sections are: (1) find experimental reaction data 
and target isotopes as accurately as possible: (y, n) data is 
preferred over (y, total) and single target isotopes are pre- 
ferred to targets with more than one isotope; (2) the cross 
section energy range should expand the energy range of in- 
terest (10-30 MeV); (3) the data is as recent as available; (4) 
data from highly relevant sources in nuclear physics was pre- 
ferred. 


D. Photon yeild 


With proper cross section data, we can apply quasi- 
absolute method by Eq. (5) if photon yield is available. Pho- 
ton yield was obtained from photon beam simulations by 
the ratio between the number of photons in a certain energy 
bin and total incident electrons. The number of photons in 
different energy bins with different sizes was counted by the 
function “hist” in statistical software R [38]. We tabulated the 
cross section data and photon yield and summed the products 
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of them to obtain the reduced reaction rate ZAR and its corre- 
sponding uncertainty. Table 1 gives an example of this pro- 
cess. In Table 1, E is the energy, o is cross section, Ag is un- 
certainty of cross section, “Energy Range” is the energy range 
centered on E, N,(£) is the number of photons in the energy 
range, Y (E) is photon yield in the energy range, Y (E)o (E) 
is the reduced reaction rate in the energy range, AY(E)o(E) 
is the uncertainty of reduced reaction rate in the energy range, 
and Ip is reduced reaction rate. After having obtained the re- 
duced reaction rate Ip and its uncertainty Alp, we inserted 
them into Eq. (6) to calculate the estimated uncertainty for the 
quasi-absolute method. 


IV. RESULTS AND DISCUSSION 


Table 2 shows the experimental results from the quasi- 
absolute method. In Table 2, h is natural abundance, E is 
the energy line we chose for calculation, T) 2 is the half-life 
of the product nuclides, Branch% is the absolute intensity of 
the energy lines, 77 is the efficiency of the detector of that en- 
ergy line, J is the reaction rate which was calculated from 
the simulation and the cross section database, P is the net 
peak area. c is the NIST certified value of concentration, 
C4 is the concentration calculated from the quasi-absolute 
method, and “accuracy” is the discrepancy of the two val- 
ues. The uncertainty of the NIST certified value is expanded 
uncertainty, with a coverage factor of k = 2 (approximately 
95% confidence) [14]. The uncertainty from quasi-absolute 
method is calculated by the standard uncertainty propagation 
formulas. 

From Table 2, one can see that the calculated value is close 
to the certified value for Ca, Mn, Co, Ni, and Zn; but apart 
for As, Rb, Sb, and Ce. The calculated concentrations are 
within a factor of two from the known concentrations. For 
some elements (e.g. As, Rb, Sb, and Ce), the results from 
the quasi-absolute method are considerably lower than the 
known concentration, close to a factor of two. This may 
stem from the fact that the cross section data we applied to 
those isotopes is an IAEA-tabulated total photonuclear cross 
section. For several elements (e.g. Ca, Mn, Co, Ni, and 
Zn), the calculated concentrations are reasonably close to the 
known concentrations (around 20%). For Zinc, two different 
reactions—(y, n) and (y, p)—reach very similar results. In 
addition, our results from the quasi-absolute method are sys- 
tematically lower, which could mean that the simulated pho- 
ton flux is higher than the real one. Many experiments have 
shown that simulated flux of standard Monte Carlo codes is 
higher than the real experimental flux by around 20%, ac- 
cording to different experimental setups and simulation pro- 
grams [39-41]. This difference may come from a combina- 
tion of several factors, such as beam emittance, uncertainty of 
beam current, beam loading, beam wandering, cross sections, 
energy dissipation in experimental setup, etc. 

Because of the limitations in measurement of experimen- 
tal parameters, computer simulation of the Bremsstrahlung 
photon beam, and cross section data, the uncertainty of the 
quasi-absolute method may be still too high for practical use 
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TABLE 1. Reduced reaction rate and its uncertainty for 7 n(y,n)°Zn reaction 


E (MeV) o (mb) A5 (mb) Energy range (MeV) N,(E) Y(E) Y(E)o(E) AY (E)o(E) 
10.82 3.5 0.1 10.82~10.92 642 0.0000642 2.247E-31 6.42E-33 
11.02 4.8 0.1 10.92~11.12 1271 0.0001271 6.1008E-31 1.271E-32 
11.22 73 0.1 11.12~11.32 1194 0.0001194 8.7162E-31 1.194E-32 
11.42 8.7 0.2 11.32~11.52 1188 0.0001188 1.03356E-30 2.376E-32 
23.82 16.1 5.1 23.72~23.92 151 0.0000151 2.4311E-31 7.701E-32 
24.02 19.8 5 23.92~24.12 147 0.0000147 2.9106E-31 7.35E-32 
24.22 13.8 5.1 24.12~24.22 68 0.0000068 9.384E-32 3.468E-32 

S Ir =1.47582E-28 + 8.08109E-30 
TABLE 2. Results from quasi-absolute method for NIST sample 

Reaction h% E (keV) Ti /2 Branch% n I(s™>) P C (ug/g) Ca (ug/g) Accuracy 

BCa(y, n) Ca 0.187 1297.09 4d12h51m40s 71.00 0.01950 5.79E-14 26660 + 290 58400 + 1900 47700 + 8300 —18.3% 

SMn(y,n) Mn 100 83485 312d6h26m40s 99.98 0.03002 5.86E-14 302104187 790 + 44 786 + 98 —0.5% 

YCo(y,n) Co 100 810.775 70d20h33m20s4 99.45 0.03089 6.51E-14 2280 +310 17.93+0.68 13.7+2.5  —23.5% 

S8Ni(y, n) Ni 68.08 1377.63 1d11h36m40s 81.7 0.01839 2.20E-14 763 + 39 81.1 +6.8 88 + 12 8.3% 

SZn(y,n)%Zn 27.90 1115.55 244d5h6m40s 50.6 0.02260 7.31E-14 19760 +290 4800 +270 3660 +420 —23.7% 

SBZn(y,p) Cu 18.75 184.577  2d13h50m 48.7 0.12344 5.24E-15 30300 + 1100 4800+270 3610+380 —24.8% 

15As(y,n) As 100 595.85 17d18h23m20s 59.40 0.04174 9.95E-14 352094206 1155439 73.747.9  —36.2% 

SRb(y,n)*Rb 72.17 881.61 32d18h23m20s 68.98 0.02846 1.44E-13 6560+230 51.041.5 28.94+3.2 —43.4% 

38 b(y,n)!?Sb 42.79 564.119 2d17h21m40s 70.67 0.04403 2.88E-13 8240 +350 45.4+1.4 249428 —45.1% 

H0Ce(y, n) Ce 88.45 165.864 137d14h46m40s 79.9 0.13253 3.16E-13 18080 +270 546422 322434 —41.0% 


in precise radiochemistry and nuclear chemistry. In addition, 
the systematically higher photon flux of simulation is a source 
of systematic error of the quasi-absolute method. 


V. CONCLUSION 


We integrated the calculation of absolute measurement 
with the simulation of the bremsstrahlung photon yield and 
raised a question whether the quasi-absolute method could be 
developed and applied to practical radioanalytical chemistry. 
The experimental results show that the quasi-absolute method 
can be effective for some elements if the cross section data is 
accurate enough and simulated photon flux is close to a real 
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